NASA 

Technical 

Paper 


April 1982 


NASA 

TP 

Vj 1979 
- . c . 1 

'■■■■■I 




Harold N. Murrow, 
William E. McCain, 
and Richard 


LOaS C\)PV. tiij 

■ ' y';'; ,v.^ 



TECH LIBRARY KAFB, NM 


NASA 

Technical 

Paper 

1979 

1982 


TECH LIBRARY KAFB, NM 



□ □b7B31. 


Power Spectral Measurements 
of Clear- Air Turbulence 
to Long Wavelengths for 
Altitudes up to 14 000 Meters 


Harold N. Murrow, 
William E. McCain, 
and Richard H. Rhyne 

Langley Research Center 
Hampton, Virginia 


IVJASA 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Branch 



r 


SUMMARY 

Measurements of three components of clear-air atmospheric turbulence have been 
made with an airplane incorporating a special instrumentation system to provide accu- 
rate data resolution to wavelengths of approximately 12 500 m (40 000 ft). Flight 
samplings covered an altitude range from approximately 500 to 14 000 m (1500 to 
46 000 ft) in various meteorological conditions. This report presents individual 
autocorrelation functions and power spectra for the three turbulence components from 
43 data runs ta)cen primarily from mountain-wave and jet-strecim encounters. The 
flight location (Eastern or Western United States), date, time, run length, intensity 
level (standard deviation), and values of statistical degrees of freedom for each run 
are provided in tabular form. 


The present data support previous data confirming that the Von Karman model is a 
realistic representation of atmospheric turbulence at wavelengths shorter than about 
1000 m ( «3000 ft). However, for wind-shear and mountain-wave situations, the spec- 
tral shape at long wavelengths (especially for horizontal components of gust 
velocity) does not conform to the Von Karman expression. It, thus, appears that 
alternate models should be considered that take into account the fact that the over- 
all spectra contain, in addition to a Von Karman component, additional power at the 
longer wavelengths. The data presented in this report, along with detailed meteoro- 
logical descriptions for each sampling flight presented in previously published 
reports, should provide adequate information for detailed meteorological correla- 
tions. Some time histories which contain predominant low-frequency wave motion are 
also presented. 


INTRODUCTION 

The required structural strength of portions of present-day transport airplanes 
is, in many cases, dictated by loads predicted to arise from atmospheric-turbulence 
encounters during the lifetime of the airplane. Presently favored turbulence design 
techniques (particularly in the United States) include the use of power spectral- 
density (PSD) methods (see, for example, ref. 1) which, in turn, require a mathemati- 
cal model to describe atmospheric turbulence in power spectral form. The model pres- 
ently used for this purpose, in most cases, is the so-called Von Karman turbulence 
model which essentially describes the power spectral density or "power" of the turbu- 
lence as a function of wavelength or frequency. (See ref. 1.) A recognized short- 
coming in the use of the Von Karman model has been the lack of experimental verifica- 
tion of the power in the long-wavelength region. The advent of new, large, flexible 
high-speed airplanes such as supersonic transports has made the long-wavelength 
region more important. Accordingly, the National Aeronautics and Space 
Administration initiated a program aimed at providing accurate power spectral 
measurements of atmospheric turbulence to wavelengths of the order of 20 000 m 
( »65 000 ft) in an attempt to describe better the power content of atmospheric 
turbulence. (See ref. 2.) 


The variable in the Von Karman model which determines the "knee frequency," or 
wavelength at which the PSD curve breaks over or flattens, is the integral scale 
value (or simply "scale") L, sometimes thought of as a maximum average eddy size. 
The value of L is the controlling factor of the model shape at long wavelength. 



The value of L = 762 in (2500 ft) , presently specified in military design require- 
ments for altitudes above 762 m (2500 ft), has been considered to be conservative and 
not solidly based on experimental evidence. In contrast to the uncertainty surround- 
ing the long-wavelength region, a considerable body of experimental evidence has 
confirmed the validity of the Von Karman model in the shorter wavelength region. In 
this region, which encompasses the short period, Dutch roll, and structural mode 
responses of subsonic airplanes, the PSD drops off in proportion to the inverse wave- 
length 1/\ raised to the -5/3 power. 

The three principal reasons why ejqjerimental data have not been available to 
date to validate the Von Karman model at long wavelengths are given as follows: 

First, rather long samples of continuous turbulence are required in order for 
the power spectral estimates to be statistically significant when data are processed 
to long wavelengths (i.e., low frequencies). Samples of sufficient duration are 
difficult to acquire. In addition, until recently it was believed that even rela- 
tively minor changes in the turbulence-intensity level during the data interval could 
not be tolerated. This concern has been alleviated by the work of reference 3. 

Second, a great deal of difficulty has always been encountered in making suf- 
ficiently accurate motion measurements of the sampling airplane at very low fre- 
quency. Low-frequency motions, arising from the pilot's actions in controlling the 
airplane and/or from the changing horizontal wind field, as a function of distance or 
time which result in an apparent drift or trend in the motion time histories - if not 
correctly measured and accounted for - can result in large errors in experimental 
power spectral values in the long- wave length region. 

Third, the lack of adequate experimental PSD measurements at long wavelengths is 
due to the nonuniform (and perhaps, in some cases, improper) use of filtering tech- 
niques between different experimenters which affect the low-frequency part of the 
time histories and, thus, the long-wavelength end of the PSD. One type of filtering 
commonly employed, referred to as "prewhitening," flattens the spectrum (i.e., makes 
it similar to a white-noise spectrum) in order to minimize "spectral window" bias 
errors (ref. 1) during a stage of data reduction; and then at a later stage the 
resulting spectrum is corrected for the effect of the prewhitening filter. It was 
demonstrated in reference 4 that such a procedure distorts the power curve in the 
long-wavelength region of the PSD and, therefore, should not be used for processing 
the data to long wavelengths (or low frequency) where very narrow spectral windows 
are used. Various methods have also been proposed and used to filter or remove from 
the time histories the trends known to have been introduced either by gyro drift or 
by the accelerometer-integration process used to determine the linear airplane 
motions. Such detrending, without extreme care in application, can filter out real 
data as well as unwanted trend data in the long-wave length range. Some early experi- 
menters believed that the removal of such trends affected only the power spectral 
estimate at zero frequency. This belief was later shown to be incorrect, however, 
because of the finite width of the spectral window, or mathematical filter, used to 
obtain the PSD. (See ref. 5.) 

The purpose of this report is to present power spectra of measured gust veloci- 
ties for a variety of meteorological conditions, with emphasis on the long-wavelength 
region, and an assessment of the applicability of the Von Karman model when this 
region is of importance. Several earlier reports are complementary to this report. 
References 4 and 5 address data-processing and measurement-system adequacy aspects. 
Two special publications (refs. 6 and 7) give some general and detailed analyses of 
data contained herein. References 8 and 9 describe the meteorological and opera- 
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tional aspects of the turbulence-sampling missions from which these data were 
obtained. Spectral characteristics related to meteorological conditions will be 
noted in this report; however, detailed statistical or meteorological analyses are 
not given. It is believed that with this report and references 8 and 9, significant 
meteorological correlations and analyses can be made according to specific require- 
ments. Power spectra are presented from 43 data runs (together with their autocorre- 
lation functions) obtained from atmospheric turbulence generated by 14 different 
specific meteorological conditions (flights). Some of the more unusual time his- 
tories are also presented (generally associated with a mountain wave or extremely 
large wind shear) in order to illustrate the nature and source of high power at the 
very long wavelengths. 


SYMBOLS 


Values are presented herein in the International System of Units (SI) and, where 
considered useful, also in U.S. Customary Units. Measurements and calculations were 
made in U.S. Customary Units. 


f 

L 

R 
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At 


u,v,w 


V 

V 
a 
B 
X 

a 

2 

a 

Q 


f requency , Hz 

integral scale value, m (ft) 
autocorrelation function 

normalized autocorrelation function, R(r)/R(0) 
lag distance, or spatial lag, km (ft) 
time-sampling interval, sec 

longitudinal, lateral, and vertical components of gust velocity, 
respectively, m/sec (ft/sec) 

true airspeed, m/sec (ft/sec) 

average true airspeed during run, m/sec (ft/sec) 

angle of attack 

angle of sideslip 

wavelength, m or km (ft) 

standard deviation (intensity level) 

variance 

power spectral density 
spatial frequency, 2n/X 
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Abbreviations : 


ACF 

autocorrelation function 

CA 

California 


dof 

degrees of freedom 

GMT 

Greenwich mean 

time 

PSD 

power spectral 

density 

U.S. 

United States 


VA 

Virginia 



AIRPLANE, INSTRUMENTATION, AND MEASUREMENT ACCURACY 

Figure 1 is a photograph of the sampling airplane. Selection of the airplane 
was based on its ruggedness, altitude- and operating-range capability, ease of 
instrumentation installation, and availability. Other details concerning the air- 
plane are given in reference 8. 

The instrumentation system and measurement accuracies are described in refer- 
ence 10. Principal airflow measurements are a, P, and V; and incremental values 
measured from their mean are used. Other measurements, which defined the airplane 
motions, were applied to obtain the corrected airflow values. Estimated accuracies 
for measured quantities are given in reference 10. Data-reduction procedures are 
given in reference 5, which also contains results from an in-flight assessment of the 
effect of integrated system accuracies on gust-velocity time histories. 


DATA ACQUISITION AND PROCESSING 

A total of 46 flights with a total of 77 data runs were made over an 18-month 
period. The flights were conducted in two geographical areas, referred to as Eastern 
and Western United States (U.S.). A total of 30 flights were conducted in the 
Eastern U.S., and 16 flights were performed in the Western U.S. For the Eastern U.S. 
flights, the sampling airplane was based at Langley Air Force Base, Virginia; and for 
the Western U.S. flights, the airplane was operated from Edwards Air Force Base, 
California. The operating areas over which the data were acquired are shown in fig- 
ures 2 and 3. Pertinent topographical features are also shown in figures 2 and 3. 

The flight numbers pertinent to data presented herein are encircled and located in 
the general area where the flight runs were made. Note in figure 2 that flight 27 
data were obtained northeast of Washington, D.C. - in the New York City vicinity. 
After "quick-look" assessments with regard to turbulence intensity and continuity 
during a run, data were processed from a total of 43 runs from 14 flights. 
Meteorological aspects of the flights are described in references 8 and 9. 

The longitudinal, lateral, and vertical components, respectively, of gust 
velocity are defined in reference 5, and measurements that are required for their 
determination are given. (As will be discussed later, an antialiasing filter was not 
required.) The autocorrelation functions and power spectra were determined from 
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time-history data digitized in 0.05-sec increments giving a Nyquist frequency 
of 10 Hz. The values of PSD are presented as a function of inverse wavelength 
(given in cycles/m); the inverse wavelength was determined from the ratio of 


2 At 
f/V 


V 


given in 


cycles/sec ^ 
m/sec , 


where V is the average true airspeed during a run. 


The autocorrelation functions (AGP) have been normalized to remove intensity 
effects by dividing the function by its variance or its value for which there is no 
time shift or lag distance. The power spectra, however, are not normalized; there- 
fore, the area under the spectral curves is equivalent to the variance a^. Thus, 
the PSD plots shift vertically according to the intensity. 


The PSD and ACF data were all processed by using 1024 time lags. The mechanics 
of the data processing thus provide an effective spectral-window width of 0.0195 Hz 
when the Blackman-Tukey algorithm and Hann window are used. (See ref. 4.) Since the 
Nyquist frequency is 10 Hz, power estimates are provided at frequency increments of 
10/1024 = 0.009766 Hz. On the ACF plots, time lag has been changed to a spatial 
lag r by multiplying the time lag by the mean true airspeed for the run. 


Before discussing the results in detail, a factor applicable to a niimber of the 
longitudinal PSD plots will be mentioned. For many of the longitudinal-component 
cases presented, it should be noted that the flattening of the power spectrum at the 
high-frequency or short-wavelength region is a result of the use of a restrictor 
provided for the pitot-static tube for high-altitude flights. The use of two dif- 
ferent restrictors for flight operations above and below an altitude of 9100 m 
(30 000 ft) to provide the proper damping for the sensitive airspeed measurement is 
discussed in reference 10. In a number of instances, the high-altitude restrictor 
was installed because of a planned mission at high altitude. However, when the high- 
altitude turbulence was not found, in-flight diversions to a lower altitude resulted 
in unavoidable use of the wrong restrictor. Consequently, the slope is incorrect at 
frequencies higher than the point of deviation from the -5/3 slope in the PSD curves 
for the longitudinal component. 


RESULTS AND DISCUSSION 

Table I presents basic meteorological and geographical information about the 
flights. Further details about each flight are given in table TV of reference 8. 

The flight- and run-number designation is retained from the actual, sampling, flight 
designation in order to allow correlation with references 8 and 9. As shown in 
table I, 3 general classifications were selected, and 14 specific conditions were 
identified. It should be pointed out that even though a predominate meteorological 
condition prevails, other meteorological influences were present to some degree. 

(See refs. 8 and 9.) The general geographic location is identified in table I as 
either VA (Virginia) or CA (California). 

Table II presents the test conditions and some of the measured statistical data 
for all the data runs identified in table I. Table II lists, similarly to table I, 
the flight numbers, dates, and data-run numbers from references 8 and 9 for ease in 
correlation. The figure numbers refer to the corresponding ACF and PSD plots pre- 
sented herein. Specific data listed for each data run include mean altitude, run 
length in distance and time, statistical degrees of freedom (dof) for individual 
power estimates (according to the method of ref. 11), and the intensity or standard- 
deviation values for the vertical, lateral, and longitudinal turbulence components. 

A goal of the study was to acquire data that would provide at least 24 statistical 
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degrees of freedom; however, from table II it is seen that this goal was not achieved 
in a number of cases. Three cases which are not considered acceptable are presented 
to demonstrate the degree of scatter in individual turbulence power estimates result- 
ing from having low statistical degrees of freedom (flight 30, run 4; flight 39, 
run 10; flight 41, run 8). 

The autocorrelation and power spectra are presented in figures 4 to 46. The 
order of presentation follows the order of table II. The format is as follows; the 
normalized autocorrelation function is presented at the top of the figure, and the 
power spectrum is presented in the lower part of the figure. In each case, for ref- 
erence, a curve is included for an analytical autocorrelation function and a power 
spectrum defined by the Von Karman expression. (See ref. 1.) Thus, for the longi- 
tudinal component. 


® (Q) 

u 


2 2L 
a — 
u n 


1 

[1 + ( 1 .339LQ)^] 


and for the transverse components. 


$ (Q) = $ (Q) 

w V 


2 L 
a — 

TC 


1 + 8/3(1.3391.0)^ 
[1 + (1 .339LQ)^] 


where Q is the spatial frequency 2%/\. An integral scale value L of 762 m 
(2500 ft), which is presently specified in the design specification MIL-A-008861A 
(see ref. 12) for continuous turbulence analysis at altitudes above 762 m (2500 ft), 
was arbitrarily selected for the reference. It is to be noted that L as used in 
the Von Karman expressions is actually the L value associated with the longitudinal 
component and is twice the L value of the transverse components. (See ref. 1, 
p. 13.) In each case, the reference curve was translated vertically to cause the 
curve at short wavelengths to coincide with the measured data. Thus, the standard- 
deviation value of the reference is different for each individual spectrum. Each 
figure is in three parts to present in order the vertical, lateral, and longitudinal 
gust-velocity data for a given flight and run. 

As stated previously, the principal purpose of this sampling program was to 
investigate the long-wavelength air motion and measure its magnitude or power rela- 
tive to the shorter wavelength turbulence such as has been measured in several 
earlier sampling programs. Figures 4 to 26 present turbulence data categorized 
according to mountain-wave or orographic phenomena. 

Figure 4 presents results for a 680-sec run made at an altitude of about 1924 m 
(6311 ft) over the Appalachian Mountains in Western Virginia. (See table II.) The 
number of statistical degrees of freedom is 27. In figure 4(a), it can be seen that 
at the higher frequency or shorter wavelengths the variation of power density with 
frequency exhibits the -5/3 slope characteristic of most previous measurements and 
corresponds to the prediction of the Von Karman ejqsression. Below an inverse wave- 
length of about 10 cycles/m, however, the data deviate significantly from the ref- 
erence Von Karman curve. In this region, the power is less than the reference curve 
by roughly a factor of 2 for an inverse-wavelength interval of about 10”^ cycles/m. 
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and then the slope or power increases and the curve crosses over the reference curve 
for the longer wavelengths measured. This increase in power at long wavelengths 
would be ej{pected from observation of the normalized autocorrelation function which 
is presented at the top portion of the figure. The low point of this autocorrelation 
function, which extends well into the negative region and occurs between a lag dis- 
tance r of 6 km ( 19 685 ft) to 7 km (22 966 ft), is typical of cases where atmo- 
spheric turbulence is superimposed upon a strong wave. The autocorrelation function 
for the longitudinal components shown in figure 4(c) has a similar appearance because 
the run was made perpendicular to the waves. The minimum occurs at one-half wave- 
length of the strong wave. This can be verified by examination of the time history 
which is presented in figure 47(a). (Some significant time histories are discussed 
in the appendix and are presented in figs. 47 to 49.) Here, it is seen that the 
average transverse time of 1 cycle of the mountain wave is about 1 V 2 niin, which pro- 
duces a one-half wavelength of approximately 6.3 km (20 669 ft) when the average true 
airspeed for the run of 139.6 m/sec (458 ft/sec) is taken into account. (An under- 
standing of the appearance of autocorrelation functions associated with mountain 
waves can be related to the fact that the autocorrelation function of a pure sine 
wave is essentially a cosine function starting at 1 at r = 0 and going to -1 at 
r = \/2, and then back to 1 at r = X. ) 

The reference curve on the ACF plot (fig. 4(c)) is for the Von Karman model 
with L =^762 m (2500 ft) . As can be seen, the experimental ACF is unlike that for 
the Von Karman model because of the predominance of the wave effect. It is to be 
noted that the one-half wavelength of 6.3 km ( =»20 600 ft) indicated by the ACF should 
translate to a strong peak in the power spectrum at about 8 x 10~^ cycles/m. How- 
ever, such a long wavelength is essentially below the resolution capability of the 
present spectral analysis. 

Several of the autocorrelation functions (those with very low statistical 
degrees of freedom) exhibit discontinuities. For example, figure 15(b) has a dis- 
continuity at r = 9 km ( «29 500 ft). For runs processed with a large number of lags 
in comparison to the run length (thus producing low statistical degrees of freedom) , 
such discontinuities are observable at a value of r equivalent to one-half the run 
length and result from "noise" introduced by the mechanics of the data-reduction 
procedure. 

Evidence of strong periodic content is discernable for all of the mountain-wave 
flights (figs. 4 to 26) from examination of the ACF in a manner similar to that pre- 
viously discussed in the example of figures 4(a) and 4(c). The strength of the wave, 
and the velocity component upon which it appears, is highly correlated with the rela- 
tive angle between the flight track and the orientation of the ridge tops. The ver- 
tical and longitudinal gust-velocity components exhibit strong-wave content for 
tracks perpendicular to the ridge tops. Little evidence of waves is present for 
tracks parallel to the ridge for any component, since waves were not being inter- 
cepted by the airplane. 

In general, the PSD of figures 4 to 26 shows that mountain-wave turbulence pro- 
duces greater long-wavelength power than the Von Karman model. This is particularly 
true for the vertical and longitudinal components when the flight path is perpendicu- 
lar to the ridge line or wave orientation. In fact, with the possible exception of 
the vertical component shown in figure 17(a), none of the mountain-wave data come 
very close to fitting the Von Karman model at low frequencies with a value of refer- 
ence L of 762 m (2500 ft) . Even for this run, the knee of the PSD does not fit 
very well. Since no strong wave is discernable on any component of this run, the 
turbulence could have been generated by rough terrain only (rather than by a mountain 
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range or ridge line) and, therefore, possibly should not be classed as "mountain-wave 
turbulence." Although the PSD of the longitudinal component for flight 34, rxin 5, 
fits the reference model quite well (fig. 19(c)), the ACF indicates a wave effect at 
X/2 of about 7 km ( «23 000 ft) which would result in power higher than the reference 
at an inverse wavelength somewhat lower than the lowest point of the PSD. 

At the higher inverse wavelengths, all of the data seem to have the -5/3 power 
drop-off of the Von Karman model, with the possible exception of the vertical com- 
ponent of flight 39 (runs 3, 5, and 7) shown in figures 22(a), 23(a), and 24(a), 
respectively. The power drop-off for these runs seems to be slightly less steep; the 
cause for this is not known. 

Another general observation is that for data which include wavelengths such as 
those presented herein, the standard deviation for the horizontal components is 
nearly always greater than that for the vertical components. This is easily observed 
from inspection of table II. Since long-scale variations are a part of the nature of 
the overall wind field which, of course, is horizontal in direction, the horizontal 
components will exhibit more overall power. This suggests that although the turbu- 
lence is isotropic at short wavelengths, as has been traditionally assumed, isotropy 
does not hold for turbulence in the atmosphere at long wavelengths. 

Figures 27 to 35 present the spectral data for wind shear including the cases 
where the jet stream is directly involved. From inspection of table 11(b), it is 
seen that the variability in standard deviation was small. The intensity of the 
vertical component was always lower than that for the horizontal components and did 
not exceed 0.9 m/sec (2.9 ft/sec), whereas the intensity of the horizontal components 
was somewhat more variable from sample to sample. 

In contrast to the mountain-wave turbulence, a number of the vertical-component 
cases contain less power in the long- wave length region than the reference curve for 
the Von Karman model. (For example, see fig. 28(a).) In every case, however, the 
horizontal components contain much greater power in the long-wavelength region than 
the reference curve. Examination of the ACF for the cases of wind shear and/or jet- 
stream effects indicates that the time histories are dominated by low-frequency drift 
or trendlike effects, which would show up in the PSD at frequencies near zero. This 

is evidenced in the ACF by a curve with low slope. The spectra all appear to follow 

the -5/3 slope of the Von Karman model in the shorter wavelength region. 

After studying the ACF of flight 38 (figs. 31(a) and 32(a)), it is apparent that 
a strong wave is present on the vertical component. The indication is, therefore, 
that the turbulence for this flight was predominately "mountain wave" in nature, 
rather than the "short-wave- trough movement" suggested in reference 9. 

Figures 36 to 46 present data that are considered to represent special meteoro- 

logical conditions that cannot be considered representative of the other categories . 

On flight 8 (fig. 36) , turbulence associated with convection was sampled at an alti- 
tude of =460 m ( 1500 ft) . This case was presented previously in reference 6 and is 

A >• 

shown to fit the Von Karman model remarkably well when a value of L of 305 m 
(1000 ft) is used for the vertical component. Somewhat larger values of L are 
required to fit the horizontal components: about 600 m ( «2000 ft) for the lateral 

component and 1200 m (=<4000 ft) for the longitudinal component. Sample time his- 
tories are shown for this low-altitude convective turbulence in reference 6. 

Although the time histories for all three components are very similar in appearance, 
the ACF and PSD indicate that even low-altitude convective turbulence is not iso- 
tropic in the long-wavelength region. The rapid drop in the power spectra at the 
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shortest wavelength end (i.e., below the -5/3 reference curve) for figure 36 is a 
result of an experiment employing a digital antialiasing filter. It was determined 
that aliasing was not a problem, and use of the filter was, therefore, discontinued 
for all other data. 

Flight 32, run 4, was originally presented as a high-altitude wind-shear case in 
reference 6. (See fig. 39.) The meteorological and other aspects of this flight 
have also been reported in detail in reference 7. The meteorological analyses indi- 
cated that the turbulence encountered on this flight was triggered by wind shear; 
however, some of the runs appeared to contain orographic effects, and the intensity 
was generally greater for those than for the remaining cases for wind shear and/or 
jet-stream effects. For that reason, flight 32 was grouped with the "special cases" 
in the present paper. The time history of run 4, which is shown in reference 6 (also 
see the appendix of the present paper) , is believed to be typical of high-altitude 
wind shear and shows no orographic or wave effects. The horizontal components are 
illustrative of the cause of the high power near zero frequency on the PSD curves, 
which is typical of wind-shear turbulence. As previously shown in reference 6, the 
vertical component for run 4 (fig. 39(a)) fits the Von Karman model quite well, indi- 
cating an L value of about 305 m (1000 ft). The horizontal components have con- 
siderable more power at the long wavelengths and do not conform to the Von Karman 
model, as is evident from the ACF curves. 


Another special case which appears to conform to the Von Karman model quite 
well, as far as the vertical component is concerned, is flight 42, run 6. (See 
fig. 44(a).) Flight 42 was conducted in mixed conditions associated with low- 
altitude thermal convective action and sea-breeze convergence. (See ref. 9.) Run 6 
was a low-altitude run exhibiting moderate and greater turbulence intensity. The 
horizontal components for these data also contain much more power at the long wave- 
lengths than the reference curve, with the ACF curves indicating that the power is at 
a very low frequency. 


A A 

Data presented herein indicate that the Von Karman model is good at short wave- 
lengths, but significant differences can be present at long wavelengths depending on 
the atmospheric situation. Some recent studies related to alternate modeling of 
atmospheric turbulence are reported in references 13 to 17. In general, the alter- 
nate models take into account the fact that, for power spectral representation to 
long wavelengths such as the measurements described herein, "turbulence" by the con- 

A A 

ventional definition may be described by a curve composed of a Von Karman component 
which is adequate at short wavelengths; but an additional term must be included to 
account for the additional power at the longer wavelengths because of the special 
long-wave content . 

The time histories presented for figures 47 to 49 are discussed in the appendix 
for cases of wave motion to point out significant characteristics that cannot be 
easily observed from power spectra. 


CONCLUDING REMARKS 

Measurements of three components of clear-air atmospheric turbulence have been 
made with an airplane incorporating a special instrximentation system to provide 
accurate data resolution to wavelengths of approximately 12 500 m (40 000 ft) . 

Flight samplings covered an altitude range from 500 to 14 000 m (1500 to 46 500 ft) 
in various meteorological conditions. Individual autocorrelation functions and power 
spectra have been presented for the three turbulence components from 43 data runs. 
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The majority of the data are from mountain-wave and jet-stream encounters. The use 
of power spectral-density design techniques for the response of airplanes to atmo- 
spheric turbulence requires a mathematical model describing atmospheric turbulence in 
power spectral form. 

The present data support previous data confirming that the Von Karman model is a 
realistic representation of atmospheric turbulence at wavelengths shorter than about 
1000 m (3281 ft), and it is thus usable for designing subsonic airplanes whose pri- 
mary response modes fall below this wavelength. (Past successful airplane design, of 
course, bears this out.) It should be noted that experimental limitations in this 
measurement program (primarily those of record or sample length) prevent the reliable 
resolution of integral scale values L greater than about 1800 m ( »6000 ft) since 
location of the "knee" of the Von Karman model for such high L values is at a 
longer wavelength than data can be reliably processed. The present data, however, 
show that although the design value of L of 762 m (2500 ft) specified in 
MIL-A-008861A, when tied to the appropriate intensity level a, results in proper 
subsonic design loads, it does not realistically represent what takes place in the 
atmosphere at wavelengths greater than 1000 m ( «3000 ft) . 

The limitations of the Von Karman turbulence model at wavelengths greater than 
1000 m ( «3000 ft) are of particular significance for many situations involving wind- 
shear and mountain-wave conditions. The difficulty in specifying an appropriate 
value of L makes its use difficult for other types of turbulence as well. 

Models which appear to be capable of realistically describing atmospheric turbu- 
lence behavior throughout the total range of wavelengths are described in NASA 
CR-145247 and CR-2913. It appears that an analytical model consisting of the sum of 
a "slow component" as a function of time, plus an amplitude-modulating function of 
time multiplied by a so-called "Von Karman component," would be required to describe 
many of the power spectra presented in this report. However, the determination of 
statistical parameters to insert into such a model to represent adequately an overall 
airplane turbulence experience would seem to be a formidable task. 

A practical approach, and possibly just as accurate, for overall design-loads 
prediction in a wavelength region greater than about 1000 m (“3000 ft) would be to 
use the Von Karman equation with a fixed and conservatively larger value of L. The 
model used in this fashion would be recognized as being approximate only, and it 
would also require additional research to determine appropriate combinations of L 
and a. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 3, 1982 
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APPENDIX 


SELECTED TIME HISTORIES 

Examples of wave phenomena can be seen in the time histories of figure 47. In 
figure 47(a), the three components of gust velocity are shown for a 14-min data run 
of a track diagonal to the ridges at an altitude of 1924 m (6311 ft) in what was 
classified in reference 13 as an Appalachian mountain-wave situation. The low- 
frequency oscillations prominent on the vertical and longitudinal time histories are 
characteristic of mountain waves. It would appear that the predominant wavelength is 
about 13 km ( »8 miles, 1.5 min). The turbulence intensity is quite variable, but it 
is interesting to note that the intensity is always high for all components on the 
positive slope of the wave. Figure 47(b) is a 9-min portion of a run on the same 
flight parallel to the ridges. Again, a significant wave pattern is obvious on the 
vertical and longitudinal components. Gust-velocity time histories from other runs 
of the same flight are shown in figures 47(c), 47(d), and 47(e). 

Time histories for a 5-min run in a mountain-wave situation at about 14 000 m 
( «46 000 ft) near the Sierras is presented in figure 48. In this case, the turbu- 
lence intensity is low, but wave motion is very predominant with a peak-to-peak 
longitudinal gust velocity of about 19 m/sec (62 ft/sec) over a 1-min interval. 
Oscillations in the vertical gust velocity were of lower amplitude and at a higher 
frequency, about 1 .5 cycles/min. 

Some other interesting wave phenomena are shown in figure 49 for a flight at an 
altitude of about 10 700 m ( «35 000 ft) that was categorized in reference 9 as being 
predominantly influenced by jet-stream effects, but in the presence of mountain 
waves. The time histories clearly show the very long wavelength (low-frequency) 
effects of shear on the horizontal components (especially noticeable in fig. 49(d)) 
with the shorter waves present on all components. The turbulence intensity ranged 
from light to very light. 
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TABLE I.- PREDOMINANT METEOROLOGICAL CONDITION AND GENERAL 


GEOGRAPHICAL LOCATION FOR FLIGHT DATA RUNS PRESENTED 


Flight 

Number of 
data runs 
presented 

Predominant meteorological 
condition 

(a) 

Flight 

location 

(a) 

Mountain-wave and/or orographic effects 

20 

5 

Appalachian Mountain wave 

VA 

28 

3 

Low-altitude wind shear and 
orographic effects 

CA 

29 

3 

Mountain-wave and low-level 
orographic effects; rotor 
zone 

CA 

30 

2 

Mountain wave 

CA 

33 

1 

Low-level orographic effects 

CA 

34 

3 

Low-level mountain wave 

CA 

39 

6 

Mountain wave and intense 
upper front 

CA 

wind-shear and/or jet-stream effects 

24 

3 

Strong vertical and horizontal 
wind shears below jet stream 

VA 

27 

1 

Low-level jet stream with 
vertical wind shear below 

VA 

38 

3 

Short-wave-trough movement 

CA 

41 

2 

Fast-moving wave trough with 
vertical and horizontal wind 
shears 

CA 

Special cases 

8 

1 

Low-altitude convection 

VA 

32 

4 

Complex orographic effects; 
vertical wind shear; jet 
stream (Death Valley flight) 

CA 

42 

6 

Low-altitude convection and 
sea-breeze convergence 

CA 


^Information for this column is taken from reference 8. 



TABLE II.- TEST CONDITIONS AND MEASURED STATISTICAL DATA 


(a) Turbulence data according to mountain- wave and/or orographic effects 


Figure 

Data 

run 

Mean 

altitude 

Data-run 

length 

Statistical dof 

®w 


<^u 

m 

(ft) 

Distance, 
km (ft) 

Time, 

sec 

(a) 

m/sec 

(ft/sec) 

m/sec 

(ft/sec) 

m/sec 

(ft/sec) 



Flight 20 on 

December 

3, 1974; 16:46 to 18:59 GMT 



4 

1 

1 924 

94.9 

680 .00 

27 

2 .22 

1 .98 

3 .79 



(6 311) 

(3 11 352 ) 



(7.28) 

(6.48) 

( 12 .42 ) 

5 

2 

1 891 

74.6 

549.95 

2 1 

1 .44 

1 .44 

2 .45 



(6 205) 

(2 44 751) 



(4.74) 

(4.73) 

(8.05) 

6 

3 

1 882 

93 .9 

686.00 

27 

1.87 

2 .34 

3 .2 1 



(6 173) 

(308 071) 



(6. 13 ) 

(7.68) 

( 10 .53 ) 

7 

4 

1 899 

55.6 

419.00 

16 

1.65 

3.59 

2 .2 4 



(6 231) 

(182 415) 



(5.42 ) 

(11.77) 

(7.35) 

8 

5 

1 920 

44.0 

333 .00 

13 

1 .59 

3 .60 

2 .90 



(6 298) 

(144 357) 



(5.22 ) 

(11.80) 

(9.50) 



Flight 28 on 

February 

14, 1975; 13:22 to 

15:31 GMT 



9 

4 

1 537 

88.0 

646.00 

25 

1 .50 

3 .07 

2 .04 



(5 042) 

(288 714) 



(4.92 ) 

( 10 .08) 

(6.68) 

10 

6 

1 897 

68.6 

496.00 

19 

1 .41 

2 .61 

2 .19 



(6 223) 

(225 066) 



(4.62 ) 

(8.56) 

(7.19) 

11 

8 

1 157 

12 4.8 

912 .00 

36 

1 .92 

4.17 

3 .94 



(3 796) 

(409 449) 



(6.29) 

(13 .67) 

( 12 .94) 

Flight 29 on February 20, 1975; 19:07 to 21:38 CMT 

12 

6 

10 432 

98.2 

537.00 

2 1 

1.10 

5.20 

3 .51 



(34 225) 

(322 178) 



(3 .61) 

( 17.07) 

(11.50) 

13 

7 

7 152 

190.8 

1018.00 

40 

1 .07 

1.92 

2 .08 



(23 464) 

(625 984) 



(3 .51) 

(6.31) 

(6.82 ) 

14 

8 

1 151 

106.3 

796.00 

3 1 

1.73 

6.43 

3.15 



(3 777) 

(348 753) 



(5.66) 

(2 1 .09) 

(10.33) 


®See reference 1 1 . 
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TABLE II.- Continued 


(a) Concluded 


Figure 

Data 

Mean 

altitude 

Data- run 

length 

Statistical dof 



<^u 

run 

m 

(ft) 

Distance, 
km (ft) 

Time, 

sec 

(a) 

m/sec 

(ft/sec) 

m/sec 

(ft/sec) 

m/sec 

(ft/sec) 



1 

5’light 30 on March 7, 

1975; 21:42 to 00 

:47 GMT 



15 

4 

14 114 

18.1 

91 .00 

4 

0 .68 

2 .56 

1.03 



(46 306) 

(59 383 ) 



(2 .23 ) 

(8.40 ) 

(3 .38) 

16 

8 

14 265 

148.7 

754.40 

29 

1.34 

5.39 

4.30 



(46 800) 

(487 861) 



(4.41) 

( 17.69) 

(14.11) 


Flight 33 on March 28, 1975; 18:47 to 20:48 GMT 


2 

2 638 

89.00 

592 .00 

23 

1 .2 7 

3 .45 


(8 656) 

(291 995) 



(4.16) 

(11.33) 


Flight 34 on April 4, 1975; 18:50 to 21:30 GMT 


3 

4 342 

55.6 

306.15 

12 

2 .05 

3 .78 


(14 245) 

(182 415) 



(6.71) 

( 12 .39) 

5 

3 874 

59.7 

332 .00 

13 

2 .34 

4.84 


(12 709) 

(195 866) 



(7.68) 

(15.87) 

7 

3 892 

88.5 

485.95 

19 

3 .82 

5.51 


(12 770) 

(290 354) 



( 12 .52 ) 

(18.09) 


Flight 39 on May 20, 1975; 18:13 to 20:26 GMT 


2 .64 
( 8 . 66 ) 

2 .80 
(9.19) 


2 

3 210 

60.6 

338.50 

13 

1 .06 

2 .3 8 


(10 532) 

(198 819) 



(3 .48) 

(7.80) 

3 

4 127 

45.4 

248.00 

10 

0 .96 

2 .48 


(13 539) 

(148 950) 



(3 .16) 

(8.15) 

5 

4 135 

86.0 

456.00 

18 

1 .08 

4.04 


(13 566) 

(2 82 152 ) 



(3 .54) 

(13 .25) 

7 

4 180 

77.7 

406.95 

16 

1 . 12 

3 .19 


(13 715) 

(2 54 92 1 ) 



(3 .67) 

( 10 .45) 

9 

4 208 

108.9 

561 .00 

22 

1 .42 

7.0 1 


(13 805) 

(3 57 2 83 ) 



(4.65) 

(23 .00) 

10 

4 157 

19.2 

100 .05 

4 

2 .08 

4.63 


(13 638) 

(62 992 ) 



(6.83 ) 

(15.19) 


4.81 

(15.79) 

4.70 

(15.43) 

4.16 

(13.65) 

6.65 
(2 1.83) 

3 .02 
(9.92 ) 

3 .65 
(11.98) 


See reference 1 1 . 


























TABLE II.- Continued 


(b) Turbulence data according to wind-shear and/or jet-stream effects 


altitude 


Data-run length 

Distance, Time, 
km (ft) sec 


Statistical dof 
for power spectra 


m/sec 

(ft/sec) 


m/sec m/sec 
(ft/sec) (ft/sec) 


Flight 24 on December 17, 1974; 16:12 to 18:04 GMT 


5 903 83.3 457.00 

(19 366) (273 294) 

5 891 98.7 535.00 

(19 328) (323 819) 

5 904 ^131.9 *^715. 00 

(432 743) 

( 19 369) '^65.1 '^353.00 

(213 583) 


Flight 27 on January 30, 1975; 15:48 to 18:07 GMT 
_ ^ 

10 605 51.8 280.00 11 0.58 

(34 793) (169 948) (1.90 

Flight 38 on April 24, 1975; 22:10 to 00:54 GMT 

10 637 63.58 337.00 13 0.75 

(34 898) (208 596) (2.47 

10 705 132.4 694.00 27 0.67 

(35 120) (434 383) (2.20 

10 649 145.7 780.95 31 0.89 

(34 938) (478 018) (2.92 


0.46 

2.31 

2.35 

(1.50) 

(7.95) 

(7.70) 

0.48 

3.04 

3.96 

(1.57) 

(9.96) 

( 12.98) 

0.71 

3.23 

3.01 

(2.32) 

( 10.59) 

(9.89) 


0.58 

1 .03 

2.23 

(1.90) 

(3.39) 

(7.30) 

GMT 

0.75 

2.83 

2.41 

(2.47) 

(9.27) 

(7.92) 

0.67 

3. 10 

1.70 

(2.20) 

( 10.18) 

(5.58) 

0.89 

5.18 

4.78 

(2.92) 

( 17.00) 

( 15.68) 


Flight 41 on May 21, 1975; 21:50 to 00:34 GMT 
34.8 184.00 7 0.' 


7 368 34.8 

(24 172) (114 173) 

7 411 46.8 

(24 314) (153 543) 


0.72 

1 .49 

1.55 

(2.37) 

(4.89) 

(5.09) 

0.79 

2.16 

2.32 

(2.60) 

(7.09) 

(7.61) 



Figure 


Data 

run 


Mean 

altitude 


Data-run length 


Statistical dof 
for power spectra 


m 

(ft) 


I Distance, 
km (ft) 


Time, 

sec 


(a) 


'^w 


'^u 

m/sec 

(ft/sec) 

m/sec 

(ft/sec) 

m/sec 

(ft/sec) 


Flight 8 on June 19, 1974; 18:09 to 19:32 GMT 


36 

2 

460 

(1 509) 

147.6 
(484 2 52 ) 

1146.3 5 

45 

1.15 
(3 .78) 

1 . 18 
(3 .86) 

1.34 

(4.41) 



Flight 32 on March 26, 1975; 19:46 to 

22 

: 19 GMT 



37 

■ 

12 989 
(42 615) 

115.8 
(3 79 92 1 ) 

634.00 

25 


2 .2 4 
(7.34) 

8.41 

(27.60) 

3 .33 
(10.92) 

38 

■ 

13 056 
(42 835) 

65.1 

(2 13 583 ) 

367.00 

14 


2 .04 
(6.68) 

4.72 

(15.50) 

3.17 

(10.40) 

39 

■ 

12 984 
(42 600) 

136.8 
(448 819) 

728.80 

28 


2 .45 
(8.05) 

7.32 
(2 4.02 ) 

4.47 

(14.67) 

40 

■ 

13 070 
(42 879} 

130 .6 
(42 8 478) 

699.00 

27 


1 .52 
(4.99) 

8.75 
(2 8 . 72 ) 

5.70 

(18.69) 



Flight 42 on June 13 

, 1975; 18:11 to 

20: 

:17 GMT 



41 

■ 

4 420 
(14 500) 

84.8 

(278 2 15) 

423 .95 

17 


1.91 

(6.28) 

2 .2 1 
(7.25) 

2 .80 
(9.20) 

42 

■ 

4 420 
(14 500) 

102 .6 
(336 614) 

512 .95 

20 


1.79 

(5.86) 

2.21 

(7.46) 

2 .07 
(6.80) 

43 

D 

4 420 
(14 500) 

57.7 

( 189 304) 

290.85 

1 1 


1.41 

(4.63) 

2 .90 
(9.53 ) 


44 

6 

1 676 
(5 500) 

79.0 

(259 186) 

544.80 

2 1 


2 .05 
(6.74) 

2 .69 
(8.81 ) 

^9 

45 

8 

3 200 
(10 500) 

54.0 

(177 165) 

2 72 . 00 

1 1 


1 .76 
(5.76) 

2 .84 
(9.33) 

1.56 

(5.12) 

46 

9 

3 200 
(10 500) 

66.5 

(2 18 176) 

334.85 

13 


2 .39 
(7.83) 

3.15 

(10.34) 

2 . 14 
(7.01) 
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Figure 2.- Location of flight tracks for flights in the Eastern United States. 
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Figure 4.- Continued. 
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Figure 7.- Concluded. 
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Figure 8.- Continued. 
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Figure 10.- Continued. 
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Figure 12.- Continued. 
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Figure 16.- Continued. 
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Figure 17.- 
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Figure 18.- 
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Figure 21.- Continued. 











•3 


i 

l/A, cycles/m 


(b) Lateral component of gust veloci 


Figure 22.- Continued 
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Figure 26.- Concluded. 
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Figure 28.- Power spectra and autocori 
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Figure 28.- Continued. 
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Figure 30.- Continued. 
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Figure 32.- Power spectra and autocorrelation functions for flight 38, run 6 
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Figure 33.- P 
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Figure 33.- Concluded. 
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Figure 34.- Continued. 
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Figure 36.- Continued. 
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Figure 37.- Continued. 
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Figure 38.- Concluded. 
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Figure 39.- Concluded. 
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Figure 41.- Concluded. 
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Figure 43.- Concluded. 
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Figure 46.- Continued. 
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(a) Flight 38, run 3. 

Figure 49.- Time histories of jet stream with long-wavelength effects at an altitude 

of about 10 700 m ( «35 000 ft). 
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